UNCLASSIFIED 


AD  NUMBER 


AD480995 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  DEC  1965. 
Other  requests  shall  be  referred  to  Office 
of  Naval  Research,  875  North  Randolph 
Street,  Arlington,  VA  22203. 


AUTHORITY 


ONR  notice  dtd  27  Jul  1971 


THIS  PAGE  IS  UNCLASSIFIED 


A  COMPARISON  OP  DUCTED  PROPELLER  THEORY 
WITH  BELL  X-22A  EXPERIMENTAL  DATA 

by 

Gary  r.  Hough  8c  Alvin  L.  Kaslcel 


TAR-TR  6510 


December  1965 


Submitted  to 

Air  Programs,  Office  of  Naval  Research 
Washington,  D.  C. 

In  partial  fulfillment  of  Contract  Nonr-4357 (00) 


Approved: 

Alfred  Ritter, 
President 


ABSTRACT 


For  the  forward  flight  regime,  a  limited  comparison  is  made  between  theo¬ 
retical  predictions  of  duct  pressure  distributions  and  data  obtained  fron.  one- 
third  and  full  scale  model  tests  of  the  Bell  X-22A  ducted  propeller  unit.  The 
theoretical  calculations  are  based  upon  studies  of  ducted  propellers  with  finite 
blade  number  which  have  been  undertaken  at  Therm  Advanced  Research,  Inc.  It  is 
found  that  the  theory  is  in  reasonable  agreement  with  experiment  and  generally 
tends  to  underestimate  the  measured  pressures.  Also,  the  characteristic  shape 
of  the  predicted  distribution  agrees  well  with  the  measured  distribution.  More 
conclusive  evaluation  of  the  reliability  and  range  of  applicability  of  the 
theory  is  needed,  not  only  for  forward  flight,  but  also  for  the  static  case  and 
in  transition. 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 
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NOMENCLATURE 


C,.t  propoller  thrust  coefficient,  TpApU^R2 

c  duct  chord,  ft 

c  duct  surface  pressure  coefficient,  (p-p  )/Hp^  ,  see  also  Kq.  (1) 

jp  00 

J  propeller  tip  advance  ratio,  U/ftRp 

p 

P  static  pressure  on  duct  surface,  lb/ft 

2 

P  free  stream  static  pressure,  lb/ft 

00 

Rp  propeller  radius,  ft 

Tp  propeller  thrust,  lb 

U  free  stream  velocity,  ft/sec 

u  total  axial  perturbation  velocity  on  duct  surface,  ft/sec 

x  axial  distance  from  iuct  leading  edge,  ft 

Ac  duct  net  loading,  or  inner  surface  pressure  coefficient  minus  outer 

p  surface  pressure  coefficient 

p  fluid  density,  lb  sec2/fti* 

ft  propeller  angular  velocity,  rad/sec 
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INTRODUCTION 


A  comprehensive  theoretical  study  of  ducted  propellers  with  finite  blade 

number  was  undertaken  at  Therm  Advanced  Research,  Inc.1-^,  the  major  portion  of 

the  study  being  concerned  with  the  determination  of  the  steady  loading  on  the 

duct  in  forward  flight  at  zero  angle  of  attack,  in  order  to  facilitate  pre- 

1  p  il 

liminary  design  applications,  the  theory  as  developed  '  '  has  been  put  into 

c  8 

the  form  of  a  simple  numerical  procedure-"  which  permits  the  rapid  calculation 
of  effects  of  varying  duct  parameters  such  as  camber,  thickness  and  chord-to- 
diameter  ratio,  as  well  as  propeller  parameters  such  as  thrust  coefficient, 
blade  number,  advance  ratio,  tip  clearance,  and  axial  location  in  the  duct. 

An  essential  part  of  this  study  would  have  been  the  experimental  evalua¬ 
tion  of  the  reliability  and  range  of  applicability  of  the  theory  as  developed. 
Unfortunately,  a  lack  of  suitable  experimental  data  on  ducted  propellers  pre¬ 
cluded  such  an  evaluation.  However,  with  the  current  design  and  development  of 
the  X-2SA  Tri-Service  VTOD  Aircraft  by  Bell  Aerosystems  Co.,  both  one-third 
scale  and  full  scale  test  results  of  a  powered  ducted-propeller  configuration 
became  available.  These  tests  were  not,  of  course,  designed  specifically  to 
provide  a  means  for  verification  of  our  theory.  However,  by  a  judicious  choice 
of  parameters  and  a  careful  examination  of  the  test  data,  it  was  possible  to 
extract  sufficient  results  to  enable  us  to  obtain  at  least  a  preliminary  com¬ 
parison  with  our  theoretical  predictions. 

The  purpose  of  this  report,  then,  is  to  make  such  a  comparison  for  the 
forward  flight  regime  at  zero  angle  of  attack.  The  first  Chapter  outlines 
briefly  the  basic  mathematical  model,  the  assumptions  made  and  the  development 
of  the  theory.  Chapter  Two  is  devoted  to  the  comparison  of  the  one-third  scale 
X-22A  test  data  with  theoretical  predictions  while  Chapter  Three  gives  corres¬ 
ponding  results  for  the  comparison  with  the  full  scale  X-22A  tests.  Finally, 
the  results  are  summarized  and  appropriate  conclusions  drawn  as  to  the  degree 
of  validity  of  the  theory. 
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CHAPTER  ONE 


BASIC  THEORY 


1.1  Formulation 

The  three-dimensional  theory  o£  ducted  propellers  in  forward  flight  as 

1-4 

developed  at  Therm  Advanced  Research,  Inc.  has  been  reported  previously 
However,  it  appears  appropriate  to  review  briefly  the  basic  ideas  underlying 
the  theory  before  proceeding  to  the  details  of  the  limited  comparison  with 
experiment.  The  mathematical  model  for  the  ducted  propeller  is  formulated  as 
follows.  A  lightly -loaded  ducted  propeller  in  uniform  motion  at  zero  angle  of 
attack  in  an  inviscid,  incompressible  fluid  is  considered.  Viewed  in  propeller- 
fixed  coordinates,  the  duct  is  replaced  by  suitable  distributions  of  sources  and 
vortices  which  are  used  to  represent  the  effects  of  thickness  and  camber.  These 
singularities  are  placed  on  an  appropriate  mean  reference  cylinder.  The  pro¬ 
peller  is  represented  by  the  classical  model  of  a  bound,  radial  vortex  line  of 
varying  circulation  for  each  blade  accompanied  by  associated  trailing  helical 
vortex  sheets.  The  path  of  these  trailing  vortices  is  determined  solely  by 
the  incoming  free  stream  and  the  angular  rotational  velocity  of  the  propeller. 

If  the  duct  geometry  is  specified  and  the  form  of  the  propeller  circulation 
distribution  is  assumed,  the  formulation  of  the  governing  equations  follows  the 
procedures  of  thin  lifting-surface  theory.  That  is,  we  require  the  flow  to  be 
tangent  everywhere  to  the  duct  surface  and  satisfy  the  Kutta  condition  at  the 
trailing  edge.  Within  the  limitations  of  the  theory,  the  ^ .  hlem  of  determining 
the  steady  duct  loading  is  equivalent  to  the  problem  of  a  ring  wing  with  the 
same  thickness  distribution  but  having  an  "ef£ectiveM  camber  composed  of  the 
original  duct  camber  plus  the  camber  required  to  cancel  the  radial  velocities 
induced  by  the  propeller  trailing  vortices  as  well  as  by  the  sources  used  to 
provide  duct  thickness. 

The  duct  source  strength  is  calculated  directly  from  the  known  rate  of 
change  of  the  duct  thickness  and  the  corresponding  induced  camber  may  be  ex¬ 
pressed  as  a  simple  matrix  operation  on  known  thickness  distribution  coeffi- 
4 

cients  .  With  regard  to  the  propeller  induced  camber,  this  is  determined  once 
the  propeller  circulation  distribution  is  specified.  The  duct  vortex  strength 
is  then  obtained  by  correcting  the  value  associated  with  the  effective  camber 
for  the  effect  of  duct  chord-to-diameter  ratio.  This  operation  has  also  been 
expressed  in  simple  matrix  form  . 

In  order  to  present  the  analytic  solution  in  a  form  suitable  for  detailed 
aerodynamic  evaluation,  certain  calculations*  requiring  the  use  of  a  digital 


*  The  propeller  circulation  was  assumed  to  be  reasonably  approximated  by  the 
Betz  optimum10  with  a  Goodman  tip  correction11  for  these  calculations. 
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computer  were  carried  out  for  representative  ducted  propeller  configurations, 
an  engineering  design  procedure  and  accompanying  charts  were  prepared  and  the 

s  Q 

results  were  combined  in  report  fomr'  .  Thus,  a  simple  procedure  is  available 
whereby  the  steady  duct  net  loading,  surface  pressure  distributions,  thrust  and 
sectional  radial  force,  pitching  moment  and  center  of  pressure  can  be  calculated 
by  hand.  Configurations  whose  parameters  fall  somewhere  in-between  or  outside 
the  tabulated  parameter  values  can  be  evaluated  by  interpolation  and/or  extra¬ 
polation  . 

1.2  Application  to  X-22A  Calculations 

The  calculation  of  the  duct  net  loading  and  surface  pressure  distributions 
was  carried  out  for  the  Bell  X-22A  one-third  scale  and  full  scale  models  at 
several  flight  conditions,  the  hand -calculation  procedures  of  Ref.  5  being  used 
for  this  purpose.  In  general,  the  X-22A  parameters  did  not  correspond  exactly 
with  those  for  which  computational  tables  had  been  prepared.  As  a  result, 
interpolation  was  required  in  several  of  the  tables  given  in  that  report.  The 
duct  thickness  and  qeoroetric  camber  distributions  were  approximated  by  seven 
and  thirteen  term  power  series  respectively. 

Two  changes  were  made  in  the  calculation  procedure  of  Ref.  5  ,  one  being  in 
the  final  calculation  of  the  surface  pressure  coefficients  and  the  other  in  the 
inclusion  of  hub  effects. 

If  u  is  the  total  axial  perturbation  velocity  on  the  assumed  reference 
cylinder  vrtiich  represents  the  duct  and  U  is  the  free  stream  velocity,  the 
full  Bernoulli  equation  for  the  surface  pressure  coefficient  cp  is  given  by 

cp  *  -2u/u  -  (u /U)2  (1) 

For  the  calculations  contained  in  this  report,  Eq.  (1)  was  used  for  cp  with 
u  composed  of  the  contributions  from  the  singularities  representing  the  duct, 
propeller  and  hub.  This  method  differs  from  that  employed  in  Ref.  5  •  There, 

either  the  well-known  linear  form  of  Bernoulli’s  equation  or  the  form  of  Eq.  (1) 

2 

with  the  term  of  order  u  containing  only  the  contribution  from  the  duct  source 
distribution  was  used. 

The  second  change  was  to  permit  the  inclusion  of  the  effect  of  the  propeller 
hub.  For  this,  the  hub  was  replaced  simply  by  an  appropriate  distribution  of 
line  sources  along  the  hub  axis.  The  effective  induced  camber  and  the  contri¬ 
bution  to  the  axial  perturbation  velocity  on  the  shroud  were  then  determined 
straightforwardly  in  accordance  with  the  methods  as  given  in,  say.  Ref.  12  .. 


CHAPTER  TWO 


COMPARISON  WITH  ONE-THIRD  SCALE  EXPERIMENTS 

2.1  Model  Data  and  Tests 

The  powered  one-third  scale  model  of  the  Bell  X-22A  ducted  propeller  was 
tested  in  the  8'  by  10'  subsonic  wind  tunnel  at  the  David  Taylor  Model  Basin. 

The  tests  were  carried  out  in  the  period  from  the  Pall  of  1963  to  the  Spring  of 
1964  and  were  designed  to  obtain  propulsion  data  in  order  to  verify  the  estimated 
performance .  Data  was  recorded  for  the  various  combinations  of  angle  of  attach, 
forward  flight  speed  and  power  input  that  were  anticipated  for  the  complete  flight 
regime  of  the  X-22A  . 

The  model  geometry  was  scaled  down  exactly,  from  its  full  scale  counterpart, 
by  a  factor  of  three  except  for  the  propeller  hub  which  was  somewhat  larger  than 
one-third  scale  in  order  to  house  the  driving  motor.  Dimensions  of  the  principal 
parameters  for  both  the  one-third  scale  and  full  scale  units  auce  listed  below. 

A  three  bladed  propeller  was  used;  further  information  regarding  the  model  and 
test  runs  may  be  found  in  Ref.  13  . 


Principal  Parameters 

1/3  -  Scale  Model 

Full  Scale  Model 

Duct  chord 

16.3” 

49.0" 

Maximum  duct  thickness 

2.8 

8.4 

Radius  of  reference  cylinder 

15.5 

46.6 

Propeller  radius 

14.0 

42.0 

Propeller  tip  clearance 

0.1 

0.4 

Propeller  axial  location 

4.7 

14.0 

Hub  length 

26.8 

43-9 

Maximum  hub  radius 

9.4 

17.5 

Static  pressure  distributions  were  measured  at  five  azimuthal  stations  on 
the  duct  with  sixteen  taps  on  the  inner  surface  and  fourteen  on  the  outer  surface. 
In  view  of  the  close  proximity  of  several  azimuthal  stations  to  the  elevons  and 
supporting  struts,  there  was  a  noticeable  variation  in  many  readings  and,  ac¬ 
cordingly,  considerable  care  had  to  be  exercised  in  selecting  what  could  be 
considered  as  reliable  data.  Consequently,  the  data  points  labelled  "experi¬ 
mental"  in  the  comparison  curves  represent  the  azimuthal  average  of  the-  readings 
retained  at  any  axial  location  on  the  shroud. 


4 


2.2  Results 

In  Fig.  1  we  see  a  comparison  of  the  theoretical  predictions  and  experi¬ 
mental  results  for  the  duct  inner  and  outer  surface  pressure  coefficients.  For 
this  case,  the  propeller  tip  advance  ratio,  J  *  U/nR^  ,  was  0.302  and  the  pro¬ 
peller  thrust  coefficient,  CT  *  T^ApU^ 7rRp  ,  was  0.129  *  values  which  corres¬ 
pond  to  customary  cruise  flight.  Despite  the  many  differences  between  the 
mathematical  model  used  in  the  theory  and  the  actual  model  used  in  the  tests, 
the  theoretical  predictions  and  experimental  measurements  are  in  reasonable 
agreement,  the  theory  generally  underestimating  the  magnitude  of  the  pressure 
coefficient.  However,  note  that  the  characteristic  shapes  of  the  two  curves 
are  very  similar. 

Variation  of  the  net  duct  loading,  Acp  ,  namely,  the  inner  surface  pressure 
coefficient  minus  the  outer  surface  pressure  coefficient,  is  shown  in  Fig.  2  . 
Here  again,  the  theoretical  curve  follows  the  experimental  data  thus  indicating 
that  the  duct  vortex  distribution  (which  is  simply  proportional  to  the  net 
loading)  is  predicted  fairly  well. 

The  calculated  pressure  coefficients  on  the  duct  and  hub  units  without 
propeller  are  shown  in  Fig.  3  along  with  the  corresponding  one-third  scale 
data.  To  illustrate  the  effect  of  the  hub,  the  contribution  of  duct  alone  is 
also  shown.  It  is  seen  that  aft  of  approximately  the  one-quarter  chord,  or 
equivalently,  the  propeller  plane,  the  hub  has  little  effect  on  pressure  co¬ 
efficient.  comparison  of  Figs.  1  and  3  reveals  that,  for  this  particular  pro¬ 
peller  thrust  coefficient,  the  propeller  contribution  is  at  most  about  10$  of 
the  total  value  of  the  pressure  coefficient. 

A  comparison  of  theory  and  experiment  for  a  larger  propeller  thrust  coef¬ 
ficient  is  shown  in  Fig.  4  for  J  =  0.242  and  cT  =  0.486  .  Reasonable  agree¬ 
ment  between  the  two  results  is  obtained,  particularly  with  respect  to  the 
characteristic  shapes  of  the  curves,  although  the  theory  still  tends  to  under¬ 
estimate  the  measured  values.  In  this  case,  the  propeller  contribution  has  now 
risen  to  about  20-30$  of  the  total  value  of  the  pressure  coefficient. 
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CHAPTER  THREE 

COMPARISON  WITH  FULL  SCALE  EXPERIMENTS 


3.1  Model  Data  and  Tests 

The  tests  of  the  powered  full  scale  model  of  the  X-22A  ducted  propeller 
unit,  conducted  at  the  NASA  Ames  Research  Center  In  the  40'  by  80'  wind  tunnel 
in  August  1964  and  March  1965/  are  reported  in  Ref.  14  .  The  overall  objective 
of  these  tests  was  essentially  the  same  as  those  undertaken  for  the  one-third 
scale  model.  In  addition,  a  model  which  had  a  modified  inlet  camber  and  duct 
diffuser  ratio  was  tested  to  determine  if  any  benefits  in  performance  could  be 
obtained  as  a  result  of  this  type  of  configuration  change.  The  essential  model 
dimensions  are  listed  in  Chapter  Two. 

For  these  tests,  the  static  pressure  distribution  was  measured  at  three 
azimuthal  stations  on  the  duct  with  eleven  pressure  taps  on  the  inner  surface 
and  eight  on  the  outer  surface.  Again,  a  selected  average  of  the  pressure 
readings  was  used  to  determine  the  experimental  points  shown  in  the  comparisons. 

Two  full  scale  test  runs  were  comp  aired  with  the  corresponding  theoretical 
predictions.  As  opposed  to  the  one-third  scale  model  test  cases,  the  elevons 
were  not  mounted  in  the  duct  unit  for  these  runs. 

3.2  Results 

The  results  for  J  «=  0.200  and  CT  *  0.276  are  shown  in  Fig.  5  and 
the  results  for  J  =  0.249  and  CT  =  0.496  ,  in  Fig.  6  .  The  values  of  Cj 
were  calculated  on  the  assumption  that  the  propeller  thrust  is  simply  the 
difference  between  the  total  thrust  and  the  shroud  thrust,  both  of  which  were 
measured. 

Generally,  the  agreement  between  theory  and  experiment  is  quite  good  for 
the  outer  surface  but  not  as  favorable  for  the  inner  surface  where  the  dif¬ 
ferences  at  several  stations  are  as  great  as  4o£.  However,  the  characteristic 
shape  of  the  theoretical  curve  again  matches  well  with  the  test  results. 

Further  calculations  at  higher  propeller  thrust  coefficients  show  that 
the  discrepancy  between  theory  and  experiment  increases  with  increasing  thrust 
coefficient.  This  is  to  be  expected  since  the  linearized  theory  begins  to 
break  down  at  large  thrust  coefficients. 

3.3  Modifications  to  Theory 

Two  modifications  to  the  theory  were  investigated  in  order  to  determine 
if  better  agreement  with  test  data  could  be  obtained.  First,  the  duct  boundary 
layer  displacement  thickness  was  estimated,  based  on  the  known  lnviscid  velocity 
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field,  and  then  added  to  the  duct  thickness  to  yield  a  new  "effective"  duct 

thickness.  The  calculation  of  c  was  then  carried  out  as  before  but  based  on 

P 

this  new  thickness  distribution.  Zn  view  of  the  extremely  small  viscous  layer, 
the  results  were  virtually  unchanged,  thus  indicating  chat  inclusion  of  the 
boundary-layer  thickness  in  the  previous  calculations  is  not  warranted. 

secondly,  we  investigated  the  effect  of  duct  reference  cylinder  location 
with  respect  to  the  propeller  tip.  That  is,  for  a  duct  of  large  thickness  — 
as  was  used  in  the  X-22A  —  it  appears  more  appropriate  to  distribute  the 
duct  vortices  on  reference  cylinders  more  nearly  approximating  the  true  inner 
and  outer  duct  surfaces.  This  was  in  fact  carried  out  for  run  #39-29  vrtiere  we 
found  that  a  relocation  of  the  reference  cylinder  produced  better  agreement 
between  theory  and  experiment  for  the  inner  surface  pressure  distribution.  This 
is  particularly  so  ahead  of  the  propeller  plane  where  the  theoretically  obtained 
values  were  increased  by  approximately  20 £  .  Aft  of  the  propeller,  the  effect 
of  relocating  the  reference  cylinder  diminishes  rapidly.  The  entire  outer  sur¬ 
face  pressure  distribution  is  relatively  unaffected  by  appropriate  relocation 
of  the  mean  reference  cylinder. 
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CONCLUSIONS 


As  a  first  step  in  evaluating  the  applicability  of  the  three-dimensional 
theory  of  ducted  propellers  which  had  been  developed  at  Therm  Advanced  Research, 
a  comparison  was  made  between  theoretical  predictions  of  the  duct  surface  pres¬ 
sure  distributions  and  appropriate  data  obtained  from  ono-third  and  full  scale 
model  tests  of  the  Bell  X-22A  aircraft  for  the  forward  flight  regime.  Despite 
the  very  limited  quantity  of  suitable  data  for  comparison  purposes,  it  was 
still  possible  to  draw  the  following  conclusions: 

The  theoretical  predictions  are  in  reasonable  agreement  with  experi¬ 
ment  and  generally  tend  to  underestimate  the  measured  pressures. 

The  Characteristic  shape  of  the  predicted  distribution  agrees  well 
with  the  measured  distribution. 

It  is  important  to  determine  accurately  the  contribution  to  the 
pressure  coefficient  of  the  duct  separately  as  this  comprised  from 
60  to  90  percent  of  the  total  value. 

Improved  agreement  between  theory  and  experiment  was  obtained  when  — 
to  make  allowances  for  the  relatively  thick  duct  —  mathematical 
singularities  were  distributed  on  individual  reference  cylinders 
representing  the  inner  and  outer  surfaces,  rather  than  on  a  single 
mean  reference  cylinder. 

In  order  to  consider  all  regions  of  practical,  interest,  more  con¬ 
clusive  evaluation  of  the  reliability  and  range  of  applicability  of 
the  theory  is  needed,  not  only  for  forward  flight,  but  also  for  the 
static  case  and  in  transition. 
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